Anisotropic scattering in angular-dependent magnetoresistance oscillations of quasi-2D and 
quasi-lD metals: beyond the relaxation-time approximation 
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The electrical resistivity for a current moving perpendicular to layers (chains) in quasi-2D (quasi- ID) metals 
under an applied magnetic field of varying orientation is studied using Boltzmann transport theory. We consider 
the simplest non-trivial quasi-2D and quasi- ID Fermi surfaces but allow for an arbitrary elastic collision inte- 
gral (i.e., a scattering probability with arbitrary dependence on momentum-transfer) and obtain an expression 
for the resistivity which generalizes that previously found using a single relaxation-time approximation. The 
dependence of the resistivity on the angle between the magnetic field and current changes depending on the 
momentum-dependence of the scattering probability. So, whereas zero-field intra-layer transport is sensitive 
only to the momentum-averaged scattering probability (the transport relaxation rate) the resistivity perpendicu- 
lar to layers measured in a tilted magnetic field provides detailed information about the momentum-dependence 
of interlayer scattering. These results help clarify the meaning of the relaxation rate determined from fits of 
angular-dependent magnetoresistance oscillations (AMRO) experimental data to theoretical expressions. Fur- 
thermore, we suggest how AMRO might be used to probe the dominant scattering mechanism. 

PACS numbers: 



I. INTRODUCTION 

Measurements of angle-dependent magnetoresistance os- 
cillations (AMRO) have emerged as a powerful probe of the 
electronic properties of low-dimensional metals. When the 
electric current is perpendicular to the metallic layers (chains) 
in quasi two(one)-dimensional metals, the measured conduc- 
tivity is found to oscillate as a function of the angle 6b be- 
tween the current and applied magnetic fiel d 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 . 
This dependence, and the weaker but observable variation 
with the field-angle fig in the plane perpendicular to the cur- 
rent, depends on the detailed shape of the Fermi surface and 
other physical parameters. Fits to data have been used to ex- 
tract experimental Fermi surfaces with a resolution compara- 
ble to the best available techniques * 1 V 2 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 
Recently, an AMRO study of overdoped cuprates revealed 
an anisotropic scattering rate with a magnitude that increases 
monotonically with temperature and scales with the supercon- 
ducting transition temperature 25,26 - 27 . Microwave conductiv- 
ity measurements under a strong field of variable direction 
provide additional insight, and are being used in Fermiol- 
ogy to augment the zero-frequency results 10,28,29 . With the 
more prominent role played by AMRO experimental tech- 
niques comes a greater need to examine each assumption un- 
derpinning the theoretical models used to interpret these data. 

Theoretical expressions for AMRO are most simply ob- 
tained within Boltzmann transport theory, which is valid when 
the quasiparticle interlayer transport is coherent. The calcu- 
lated magnetoresistance in strong magnetic fields is sensitive 
to the shape of the Fermi surface since electrons complete 
many cyclotron orbits before being scattered. Using AMRO 
data taken over a range of Og, <f>B and field strength, it is pos- 
sible to extract multiple hopping parameters and thus build a 
detailed topography of the intralayer Fermi surface. For over- 
doped cuprates, the results have shown good agreement with 
ARPES measurements and other determinations ' ' . 

In contrast to the close consideration given to the effects of 



complex Fermi surface shape on AMRO, electron-scattering 
effects have usually been treated simply: a single scattering 
rate, used as an additional fitting parameter, is assumed. In 
other words, the Boltzmann theory of AMRO is treated within 
the relaxation-time approximation 30 . 

Since the net change due to scattering of the electron dis- 
tribution at a given momentum depends on the distribution 
elsewhere, a general distribution function must be found self- 
consistently from the Boltzmann transport equation with the 
full collision functional include d 31 ' 32 . Within the relaxation- 
time approximation, the collision functional is assumed to be 
proportional to the nonequilibrium part of the electron distri- 
bution, with a proportionality constant defined as the current 
relaxation rate, or inverse current lifetime. The distribution is 
then easily found from the simplified transport equation. (The 
relaxation-time approximation is equivalent to neglecting ver- 
tex corrections in a Kubo calculation of the conductivity 33 .) 
The electron distribution obtained this way is rarely a solu- 
tion to the full Boltzmann equation, so the calculated trans- 
port coefficients are not always reliable. Even if experimental 
properties can be explained using this approach, the physical 
meaning of the lifetime extracted from fits to data is not al- 
ways clear. 

It is well known that lifetimes of electrons in a given mate- 
rial determined from different experiments on the same sam- 
ple can vary widely. For example the scattering rate deter- 
mined from zero-field transport, the transport relaxation rate 
t~} , is expected to be smaller than that determined from spec- 
tral or quantum oscillation measurements, the quasiparticle 
scattering rate r , because small-angle scattering processes 
contribute only to the latter 40,41 . This difference persists to the 
lowest temperatures if inhomogeneity with a large character- 
istic length scale is present in the sample, as may be the case 
for many low-dimensional metals 34,35 . 

AMRO is one measurement where the physical interpreta- 
tion of an extracted lifetime might not be immediately clear. 
For, although it is a transport measurement, the fact that elec- 



trons are bound in tight electron orbits and can be transferred 
from one orbit to another via small-angle scattering suggests 
that scattering processes which play no role in zero-field trans- 
port might be significant. Notably, the precise meaning of life- 
times extracted from quantum oscillations in magnetotrans- 
port and from cyclotron resonance experiments is still a mat- 
ter of discussion in the literature 28 ,36 i 37 i 38 i 39 i 40 i 41 If a lifetime 
determined from AMRO data using a relaxation-time approx- 
imation is to be usefully compared with those obtained else- 
where then its physical significance must be clarified 11,26 ' 36 . 
Addressing this issue is the main purpose of this paper. 

We calculate AMRO for simple model (quasi ID and 2D) 
Fermi liquid metals using Boltzmann theory with the full elas- 
tic collision functional included. Our purpose is (i) to eluci- 
date the meaning of lifetimes extracted from AMRO by trac- 
ing their origin from the collision functional, and (ii) to deter- 
mine whether new information can be extracted from AMRO 
by going beyond the relaxation-time approximation. It turns 
out that a single effective scattering rate does not fully account 
for the effect of scattering on AMRO but rather, in a crude 
sense, the effective scattering rate relevant to AMRO depends 
on the angle 6b between the field and current. Also, we find 
that AMRO can provide a description of the in-layer momen- 
tum (and corresponding spatial) dependence of the scattering 
cross-section that is not available from other techniques. Thus 
AMRO may prove useful for studying an unknown scattering 
mechanism. 

We here briefly summarize our analysis and results before 
providing details. The collision functional for a quasi-2D 
metal that is isotropic in the (x-y) plane is characterized with 
real parameters A„, where n is an integer, defined by 



A„ — A c 



■I 



dq z d<pP(\q z \,\<p\)d ~ e^"* 



(1) 



where the integral is over the (quasi-cylindrical) Fermi sur- 
face and P(\q z \, \<f>\) is the probability per unit area for scat- 
tering an electron between points on the Fermi surface sepa- 
rated by q z ,(f>. The scattering probability depends only on the 
magnitude of the change of each cylindrical momentum com- 
ponent. Each A n , having the dimensions of inverse time, is a 
scattering rate of potential physical significance. The trans- 
port relaxation rate is t~) = Aq whereas the total quasiparti- 
cle scattering rate is t _1 = Aco. The A„ for finite n contain 
additional information about the dominant scattering mech- 
anism. For example, the scale n over which A„ approaches 
Ace is roughly equal to the largest length scale over which the 
scattering potential varies within an x - y layer, measured in 
lattice constants (see Appendix). This type of parametrization 
of the collision integral is well known, and has been used, for 
example, to describe the effect of small-angle impurity scat- 
tering on the Weiss oscillations in a two-dimensional electron 
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gas. 

Our main finding is that AMRO is sensitive to A$ when the 
field is parallel to the current, to /!«, when the field is perpen- 
dicular to the current, and to intermediate A„ at intermediate 
angles. There are two immediate implications of this result. 
First, it clarifies the meaning of the scattering rate extracted 
from AMRO fits. Second, it demonstrates that AMRO mea- 



surements provide information about the momentum-transfer 
dependence of the scattering probability that may not be oth- 
erwise available. For, as long as the scattering potential ex- 
tends over more than one layer, so that Aq + A x , it should 
in principle be possible to extract more than one A n using 
AMRO (methods for isolating individual A n are described be- 
low). This would enable a partial reconstruction of the scat- 
tering cross-section via Eq. Q] and could reveal, for example, 
whether an unknown scatterer had atomic-length scale corre- 
lations in the layer or whether it was correlated over much 
larger distances. 

In the next section we give the Boltzmann derivation of the 
interlayer conductivity for the quasi-2D and quasi-lD metals. 
In Section III we give a simple physical picture of this result. 
In the subsections of Section IV we analyze the conductivity 
in different limits, describe how individual A n could in princi- 
ple be extracted from AMRO data, and illustrate the expected 
behaviour. In the Appendix we consider the simple example 
of random impurity scattering to illustrate the significance of 
the n dependence of A n . 



II. BOLTZMANN THEORY OF INTERLAYER 

MAGNETORESISTANCE WITH ELASTIC, ANISOTROPIC 

SCATTERING 

To calculate the resistance in a strong applied magnetic field 
we adopt the following form of the Boltzmann equatio n 31 ' 44 , 
which is lowest order in the strength of the electric field E but 
valid for an arbitrary-strength static magnetic field B, 



-TjT - !igk] = -<?E • v k 



(2) 



where g k is the non-equilibrium part of the full electron dis- 
tribution function f k , 



/k = /o(<*) + 



<9/0(6k) 



de k 



gk, 



/o(ek) is the Fermi distribution, and Vk = (de k /dk.). The colli- 
sion integral I[g k ] is of the form 



/fek] 



/ 



dS k ,P(k,k')(g k ,-g k ). 



(3) 



where P(k, k') is the probability per unit time for an elec- 
tron to be scattered from k to k', where both lie on the 
Fermi surface. The integral is over the k' Fermi surface with 
dS k > - d 2 k'/\\ k '\. This collision integral is appropriate for 
elastic scattering, which will henceforth be assumed. Well- 
known approximate extensions to a given inelastic scatter- 
ing mechanism could be used to predict, for example, dif- 
ferent temperature dependencies for different A n (as occurs 
for scattering by acoustic thermal phonons for T « To 
where Aq « (T/Td) 2 Aoo)2*. For electrical conductivity (as 
opposed to thermal conductivity) the elastic collision integral 
may be expected to give qualitatively correct results with the 
A„ temperature-dependent. 

Under the applied magnetic field electrons move along cy- 
clotron orbits on the Fermi surface (to lowest order in |E| the 



electric field can be ignored in this orbital motion) such that 
k = k(f) varies with the time variable t according to 



dk 

— = -e\ k x B, 

at 



(4) 



with h - 1 . 

To calculate the magnetoconductivity, we solve Eq. [T6lto 
obtain the momentum k(f) for an arbitrary initial value k(0) 
and insert the result into Eq. [2] to determine the distribution 
function g k (a single value of t is used for the momenta k(f) 
and k'(f) in the collision integral). Both Vk and g k depend on 
the initial momentum k(0) and vary periodically with time. 
The time-dependent current is 



1(f) 



— f 

2n) 2 J 



(2. 



dS kgk (t)\ k (t), 



(5) 



where the integral over the Fermi surface is done using the 
initial momentum k(0) coordinates. The frequency-dependent 
conductivity cr(co) (for both current and electric field along z) 
is finally obtained from the time-frequency Fourier series of 
the z component of Eq. [5] 

Below we follow this method to obtain the interlayer con- 
ductivity for the quasi-2D system (Eq. [To*b and quasi- ID sys- 
tem (Eq. l25l ). These two Eqs. are the main results of this 
article. 



A. Quasi-2D metal 

We consider a quasi-2D Fermi surface that is isotropic in 
the layer (the k x — k y or a - b plane) and weakly corrugated 
in the direction perpendicular to the layers (k z or c) with band 
energy 



£k = 



1 
2m* 



(kl + kl 



' */) ' 



2t c cos(k z c) 



(6) 



where m* is the effective mass and the last term is obtained 
from nearest-interlayer-neighbor hopping with coefficient t c 
and interlayer distance c. It is assumed that k f /2m 



» t r 



where kf is the average radius of the near-cylindrical Fermi 
surface. 

We substitute Eq. [6] into Eq. [4] with the magnetic field 
B = B(sin 8, 0, cos 6). To lowest order in v z , the c axis disper- 
sion of the Fermi surface can be ignored in the determination 
of the cyclotron trajectories and the collision integral so that 
k = (kf cos <p, kj sin <p, k z ). The z-component of the equation 
of motion Eq. @]is 



dk z 

-— = (x>ckf tan 9b sin (/> 



(7) 



and, dropping an interlayer hopping term under the assump- 
tion (kf/m*)sin((p) » t c c tan##, the intralayer components 
become 



d<p 

Tt= Uc 



(8) 



where the cyclotron frequency is 

wc = (eB/m*) cos Ob- (9) 

Eqs. [7] and [8] are solved to give 

v z (t) = 2t c c sin k z c - kfC tan 6b[cos((/) + Lj c t) - cos (/>] , (10) 

which can be used in the Boltzmann equation, Eq. [2] The 
omitted interlayer hopping term is typically small, significant 
only for electrons moving in the direction of the magnetic field 
when the field is nearly parallel to the layers. This term gives 
rise to a small coherence peak in the resistivity at Ob - n/2 in 
strong magnetic field but otherwise has little effecti ^ 10 ' 44 

The distribution function and velocity are expanded in a 
Fourier series over both momentum variables at t = 0: 



8k 



z 



gmn(t)e 



imk,(0)c+in(f}(0) 



(ID 



,(k) = ^ u„ m (t)e 



imk,(0)+imp(P) 



To henceforth avoid such cluttered notation, we use the sym- 
bols k z and <p to refer to the initial values k z {0) and <p(0) and 
indicate t dependence explicitly. 

In this simple isotropic model, the scattering probability 
P(k', k) depends only on the change in intralayer momentum 
\<p' - 4>\. Also, to lowest order in v z the scattering probability 
may be treated as a function of \k' z — k z \. The Fourier expansion 
thus diagonalizes the collision integral, which can be written 



I[gk\ = - 2^ gmn(t)A„ 



^imk-c+in<p 



where the collision parameters A„ m , defined by 



Xn/c pi 

dq z 
n/c JO 



d4>P{\q z \, |fl)(l - e 



_ imq-c+itup 



(12) 



(13) 



have the dimensions of inverse time. 
The Boltzmann equation Eq. |2]becomes 



(t- + Amnjgrrm = -eu mn (t)E(t), (14) 



which has the solution 



n(t) 






dt'e 



imn(f-t) 



u mn (t')E(t'). 



(15) 



The m Fourier number corresponds to a label of the layers 
in real space, i.e. the set of coefficients g m <„, for all n and a 
particular m', describes the difference in the electron distri- 
bution between layers separated by m! lattice constants along 
the c crystal axis. Since the band energy Eq. |6]contained only 
nearest-plane hopping terms (and interlayer motion within a 
cyclotron orbit was excluded by dropping terms of order v z ) 
only u\ n and u-\„ terms are nonzero. So, only the m — ±1 



terms in the distribution function are present and we can de- 
fine single-subscript collision parameters A„ = A ln = A-i„, as 
first introduced in Eq. Q] 

To carry out the <p and t integrals of the nested sine func- 
tions appearing in Eqs. [14] and [3] it is convenient to intro- 
duce the identity e '' sin(v) = E^-« J m (z)e ix , where J m (z) is an 
rath order Bessel function of the first kind, twice into the both 
the expression for gk(0 and v z (t). (Note, in the quasi-2D and 
quasi- ID models, the relationship between the corresponding 
components of the conductivity and resistivity tensors is sim- 
ply Pzz - io'zzf i we henceforth denote cr zz by cr ± .) The final 
result for the conductivity is 



o-_l(A>) 
0"j_O^O 



= Z ^W^WrTTr 

„ „__™ A n 1 + I — 



, ui-poic \2 



where the argument of the Bessel function 

y - kfCtanOs, 
the cyclotron frequency oj c is 

a>c — too cos 6b, 



(16) 



(17) 



(18) 



where a>o - (eB/m*) and cr±o is the zero field c-axis D.C. 
conductivity, given by 



0\lo = 



2e l t l c ck f 1 
nvf Aq 



(19) 



B. Quasi ID metal 



where v - s,d and u s mn 
v+(k) + v:(k) and vj(k) 
parameters are given by 



and u mn are the Fourier series for 
- v"(k), respectively. The collision 



I 



dq y dq z P(\q,\,\q z \)(l 



^imq z c+mq y t 



and 



^mn + 2 



I dq y 



dq z P,(\q z \,\q\)t 



| \ imq z c+inq y b 



(23) 



(24) 



where f(k) = ^o(k) + Pi(k) is the total scattering probabil- 
ity per unit time, Po(k) the intrasheet and Pi(k) the intersheet 
component. Just like for the quasi 2D case, only ra = ± 1 terms 
are important so we define single-subscript collision parame- 
ters XL = A v . = i v , . 

The remainder of the calculation is identical to the 2D case. 
The final result is 



o~±oAq 






»7v 



+ ( 



Ci>—pci)c \2 



(25) 



where t] s - 1 + (-1)" and t/^ = 1 - (-1)". The argument of the 
Bessel function is y = 2(tf,c/vf) tan 0^ and the prefactor 



c\i_o/lo - 



2e 2 f c 2 c 
v/for ' 



(26) 



Note that the presence of the factor tbc/v/, which will be 
smaller than one for quasi- ID systems, in the argument of 
the Bessel function suggests that observable AMRO will be 
restricted to higher values of 9b than in the quasi-2D case. 

If intersheet scattering is ignored, so A^ - A'l, then Eq. [25] 
reduces to the form of Eq. [16] 



To model a quasi- ID system we use the band energy 
e k = Vf(\k x \ - kf) - 2tb cos(kyb) - 2t c cos(k z c) 



(20) 



where both interchain hopping parameters fj, and t c are small 
compared to V/kf, The Fermi surface consists of two sheets, 
located near k x = ±kf, that are weakly dependent on k y and 
k z . Solving the Eqs. of motion, with a magnetic field B = 
B(sin 9, 0, cos 9), we obtain for the c-axis velocity on the k x — 
+kf sheets 



v*(f) = 2? c csin k z c+kfc tan 9[cos(kyb±Q)ct)-cosk y b] 



(21) 



The c-axis conductivity is calculated in much the same way 
as for the quasi 2D system. A slight complication results from 
the need to consider both intrasheet and intersheet scattering 
processes. If we write g^ for the distribution on the k x - ±kf 
sheet then intersheet scattering in the collision integral cou- 
ples g^ with g k . However, Fourier expansions of the sum 
Sk ~ §t + #k an d difference g^ — g^ - g^ diagonalize their 
respective collision integrals. The Fourier components of the 
distributions are determined from 



I n f inn J&mt 



-eu v mn (t)E(t\ 



(22) 



C. Relaxation-time approximation 

If the collision integral is treated using the relaxation time 
approximation then the n dependence of A n is ignored and a 
single transport relaxation rate, say t^ 1 is assumed. If we set 
all A n equal to r^ 1 in Eq. [l6]or[25](in the latter, no distinction 
is made between intrasheet and intersheet scattering) then, us- 
ing the Bessel function identity £m=-oo 4(x) = 1, both Eqs. 
reduce to 



o~i(") 
C\l0 



= Z J >\ 



1 



+ (b) - pCO C ) 2 T B 



2JI- 



(27) 



which has been found previously by Moses and McKenzie 7 . 



D. Contrast with the intralayer conductivity 

The intralayer conductivity component cr\\ = <j xx may eas- 
ily be calculated within the same simple models. Taking the 
electric field to lie purely along the x-axis, and repeating the 
above calculation for the quasi-2D metal, we find: 



h)\ 



o"||o^oi A 2 m +{u>- cx) C ) 2 



(28) 



where crpAoi = 4ne 2 Ekf/(mc). In this expression, there are 
no angular-dependent oscillations and the only relevant colli- 
sion parameter is the intralayer transport relaxation rate Aq\. 
All the novel effects in the inter-layer transport, Eq. [16] arise 
from the non-trivial /^-dispersion, and are absent from the in- 
tralayer transport coefficients in this model. 

For the quasi ID system, the yy component of the con- 
ductivity is of the same form as Eq. [25] (with the b and c 
axes interchanged). For the xx component, since the velocity 
v^ = ±kf/m*, the right side of Eq. [22]is equal to for v - s 



and to -leEkf/m* for v - d. 
conductivity cr\u\\{oS) is thus 

CT||(W) 



The result for the quasi- ID 



*00 



cnMoo 



«,) 2 W 



(29) 



where the prefactor is (T||o/Iq = e 2 (kf/m*) 2 /(jrbc). 

The collision parameters that enter the interlayer conductiv- 
ity A„ = A\„ describe the relaxation of a difference in current 
density on adjacent layers. They do not include any of the 
Ao„ parameters, which describe the relaxation of current den- 
sity variation within a single layer and of which the intralayer 
transport relaxation rate Aq\ is one member. Nevertheless the 
fact that Ai x = -looo - ^moo implies that the interlayer conduc- 
tivity is sensitive to the total quasiparticle relaxation rate in a 
tilted field (indeed, as discussed below, the total quasiparticle 
relaxation rate is the relevant quantity when the magnetic field 
is parallel to the layers). 




FIG. 1 : (Color online.) Interlayer current relaxation of a quasi 2D 
metal in a tilted magnetic field. The Fermi surface is shown as the 
dashed cylinder and the electron distribution in an interlayer (i.e. || 
k-) electric field as the hatched region. Electrons undergo cyclotron 
motion in the magnetic field, tilted by 6 B from the current direction. 
The two upper (blue) bundles of electrons are initially at the same 
k z and have the same density. For 8 B t their orbits separate the 
bundles in k z , which allows scattering _L k- to relax the current of 
each. Bundles on the right (blue and green), initially displaced in 
k- with different densities, orbit in phase so scattering || k- always 
exchanges electrons between them. 



III. PHYSICAL PICTURE OF SCATTERING AND AMRO 

Before proceeeding to analyze Eq. [16] in more detail, it 
is worth discussing its qualitative features (only the quasi-2D 
system for u> — will be considered). The oscillatory be- 
haviour of the magnetoresistance, which is captured by Eq. 
[27l has been discussed elsewhere^. Here we focus on the ef- 
fect of scattering and, in particular, on the difference between 
the result obtained using the full collision functional Eq. [16] 
and that found using the relaxation time approximation Eq. 
[27] The difference results from the appearance of multiple 
collision parameters A„ in Eq. Q~6] For this discussion it is 
helpful to recall that small-angle scattering within a k x - k y 
plane contributes to A„ for large n but has no effect on Aq (i.e. 
the integrand in Eq.[T]vanishes when n — and q z =0). 

The factor J 2 (y) in Eq. [16] originates from the 
Fourier expansion in (f> of the interlayer current distribution 
g(k z , <p, t)v z (k z , <p, t) at the zero of time while Jl n (y) comes 
from the current distribution at a different time. So Eq. [16] 
has the form of a current-current correlation function expected 
from the Kubo formula 45 . The "incoming" current (that asso- 
ciated with p) and "outgoing" current (associated with p — n) 
are coupled by a non-trivial vertex involving collision param- 
eters A„ . In the relaxation time approximation, the non-trivial 
vertex factors are ignored so the currents decouple and the ex- 
pression reduces to Eq. [27] 

A Bessel function J p (x) becomes small once its order p ex- 
ceeds its argument x, so the terms p that contribute to Eq. 



[16] are those for which p g kfCtanOg- This means that the 
incoming current distribution varies more rapidly around the 
cylindrical Fermi surface with increasing Ob- Since \p - n\ ^ 
kfCtanO, the outgoing current distribution also varies more 
rapidly, which enables collision parameters A„ for larger n to 
become involved. Thus AMRO becomes more sensitive to 
small-angle scattering in the k x - k y plane as Ob is increased. 
This behavior can be better appreciated using the following 
simple picture. 

In the absence of a magnetic field the conductivity is Eq. [19] 
and the nonequilibrium part of the electron distribution is 



g k = eEv z A 



(30) 



which is shown as the hatched region in Fig. [T] It is indepen- 
dent of <p, so a Fourier expansion of the current distribution in 
<p has only the p - term. Scattering perpendicular to k z con- 
nects points of equal current distribution, with no net effect, 
so only Aq (no other An) appears. 

Now suppose we start from the B - distribution Eq. [30] 
and, at t = 0, turn on a strong magnetic field at an angle 6b 
from the interlayer electric field. Electrons make cyclotron or- 
bits as shown in Fig.Q] Two bundles of electrons that initially 
had the same k z (and thus the same density gk) but different 
4> will move apart in k z as they follow their respective orbits. 
In the same way, bundles of unequal densities are brought to- 
gether onto a given k z - const plane. The current distribution 
at a given k z thus becomes dependent on at times t > 0. If 
0b is large so the orbits are elongated, this variation with 



is rapid because bundles of widely different density are being 
brought together. This is the reason why increasingly large p 
contribute in Eq. Q~6]as 6b is increased. 

If the current distribution varies rapidly as a function of <p 
then small-angle scattering within a k z — const plane can ef- 
fectively relax the current. Inspecting Eq. [TBJwe see that this 
is what happens near 6b = n/2. For, as 6b increases, terms 
in the sum for which p and p - n differ substantially appear. 
The incoming and outgoing current distributions then differ at 
points narrowly separated in <p, which allows A n for large n, 
i.e. small-angle scattering along a k x - k y plane, to affect the 
conductivity. 

While the magnetoresistance becomes more sensitive to 
small angle scattering perpendicular to k c as 6b increases, 
there is no change associated with scattering parallel to k z . 
This may be understood by following the evolution of two 
bundles of electrons initially at the same <p but different k z . 
When we turn on the field at t — 0, the bundles follow orbits 
displaced in k z while their phases </> remain equal. Scattering 
along k z transfers an electron from one bundle to the other 
just as it did at t — 0. This is true regardless of the value of 6 B . 
So scattering with momentum transfer parallel to k z is equally 
effective at all 6b- (Of course, this qualitative anisotropy in 
the effect of scattering ± and || to k z is a consequence of our 
dropping higher order terms in v, for the quasi-2D system.) 

These simple considerations are sufficient to understand the 
qualitative effect of scattering in AMRO for the quasi-2D sys- 
tem. The interlayer magnetoresistance becomes increasingly 
sensitive to small-angle scattering in the momentum plane 
perpendicular to the current as the field angle 6b is increased. 
This is why AMRO might prove a useful probe of the in- 
tralayer properties of scatterers. This physics is missed en- 
tirely in the relaxation time approximation, which treats all 
scattering processes as equivalent. 



IV. ANALYSIS OF AMRO AND THE EXTRACTION OF 
COLLISION PARAMETERS 



additional collision parameters in order to extract information 
about the scatterer. 

For the remainder of this section we analyze the behaviour 
of Eqs. [T6landl25l We are mainly interested in how A„ might 
be extracted from experiment in various limiting cases. 



A. Field along the layers 

If the magnetic field is directed along the c-axis then only 
the zeroeth order Bessel functions contribute to the sum and 
Eqs. n~6l and |25l reduce to the zero-field conductivity <x(0) = 
<x , i.e. it is sensitive to Ao alone. As 6b approaches n/2 a large 
number of terms contribute and the expression is dominated 
by large n and p. Here, A n can be replaced by A m and the 
n sum easily done. Both the pth order Bessel function and 
the Lorentzian are capable of cutting off the sum over p and 
which one actually imposes the cutoff depends on the size of 
Qc^co wriere ^c - kfc(eB/m*). The n/2 limit of Eq. [T6lis, 
for (l) « Q.c, /loo 



o- ± (cj = 0) 1 



1 



°"W() Ace ^1 + (ilc/Aoo) 2 



9 B = n/2 (31) 



(This result has previously been obtained within the relaxation 
time approximation by Schofield 44 ). The ID result is of the 
same form with D.c = 2(eB/m*)tb/(nvf) since, as seen from 
Eqs.|23]and|24Ui,=4,. 

In the Q.clAaa » 1 limit of Eq. [31] the conductivity is 
inversely proportional to field and independent of scattering, 
i.e of /loo. In the opposite limit, Qc/^oo << 1 the conductivity 
is then proportional to 1 //loo, and independent of the field. 

Comparing the 6b — and 6b - n/2 limits of Eq. [16] we 
see that the c-axis magnetoresistance is sensitive to either the 
transport relaxation rate Aq or the quasiparticle scattering rate 
/loo depending on whether the field-angle 6b is perpendicular 
or parallel to the metallic layers. Information about A n for 
finite n is available at intermediate angles. 



The accuracy of the relaxation time approximation in 
AMRO will depend on the nature of the dominant scattering 
mechanism. For scattering by atomic-scale defects, Ao and 
/loo will differ by at most a factor of order unity. If the dom- 
inant scatterer is of electronic origin then, for the anisotropic 
systems under consideration, the scattering potential will be 
expected to have only short-range interlayer correlations and, 
once again, Aq » /loo- Only when the scattering is due to long- 
wavelength phonons or other inhomogeneities extending over 
many layers will there be strong n-dependence in the collision 
parameters A„. 

Thus the relaxation-time approximation should give qual- 
itatively correct results for interlayer transport in low- 
dimensional systems for most cases of interest. But given the 
quantitative accuracy with which the Fermi surface and other 
properties have been measured using AMRO, a more precise 
description of the effect of scattering appears prudent. Also, 
as long as there is an observable difference between Aq and /l M , 
we can take advantage of the sensitivity of Eqs. [TBI and [25] to 



B. Strong-field limit (oj c » A„) 

AMRO is seen when a>c » Ao at small 6b- In this case, the 
sum over p in Eq. [TBJwill be dominated by the p — term ex- 
cept very close to the conductivity minima (where the zeroeth 
order Bessel function vanishes) and close to n/2. Everywhere 
else, the conductivity will be given by 



o-Aio = 0) %y) 



o~±oAo 



/lo 



(4(7) 
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4(7) 
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(32) 



The corresponding expression for the ID metal is 

o-j.(o> = 0) J lW>i T 2, s 4(r) 4.(7) 



(4(7) 



o~±qAq Aq \ Af A s 2 



...-), (33) 



with the A% parameters associated with the even-order and A d n 
parameters with the odd-order Bessel functions. 



C. Weak field limit (oj c « A„) 

In the opposite limit, in which ate « A„ for all n, the 9b 
dependence of the conductivity is much weaker. The second 
term in the denominator will be negligible for all p and Eq. [T6l 
becomes 



<r ± (co = 0) 

O"j_0/to 



24(r)4(rM; 



1 
p+n' 



(34) 



n,p 



This is a weighted average of A' 1 for n < n max « A;^ctan(0 B ). 
It suggests that AMRO fits using a single scattering parameter 
might be improved somewhat by allowing the scattering rate 
to depend on field-angle 9g, although this is clearly a crude 
treatment. 

For the quasi- ID case, the expression in this limit is very 
similar 



0\l(W = 0) 1 



o~±oA) 



= ? Z W(r)(W' + AK+nT 1 )- 



(35) 



D. Finite frequency conductivity 



The w-dependent conductivity provides additional routes to 
obtain the collision parameters. For a given frequency at, the 
conductivity will have peaks as a function of magnetic field 
strength when at = p'at c f° r some integer p'. So if the field 
strength B is close to B w /(p' cos 9b) where (eB^/m*) = at, the 
conductivity will be 



crjat) fyv) 



o~±oAo p' cos ( 



Y J J 2 P ,- n (y)L P ,AB) + ... 



(36) 



where the additional terms are smoothly field-dependent and 



Lp'n(B) - 



(A n lp' cos 9b) 
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p'm" cos 8b 



(37) 



If 0b is reasonably small then the zeroeth order Bessel func- 
tion in the sum over n will dominate, so n - p' and the 
field dependence is a Lorentzian with a peak located at B — 
BojKp' cos 9 b ) having width equal to A p ,/(p' cos 9 B ). Thus 
field scans of the conductivity measured at finite frequency 
and reasonably low 9b could be used to extract individual A„ 
for small n. At higher 9b the expression would no longer be 
Lorentzian, but rather a weighted sum of Lorentzians with dif- 
ferent widths A„ (n within roughly k/C tan 9b of p' contribute 
to this weighted sum). 

For the ID metal this discussion also applies with the only 
difference that the width of the corresponding Lorentzian is 
either A s /(p' cos 9 B ) or A d /(p' cos 9b) depending on whether 
p' is odd or even, respectively. 



E. Example: Gaussian scattering probability 

To illustrate the behaviour of Eqs. [16] and Eq. [25] with a 
simple example we consider a Gaussian scattering probabil- 




FIG. 2: (Color online.) Collision parameters A„ for scattering that 
favors small-angle scattering. Plotted is the n dependence of the col- 
lision parameters A„ for the Gaussian scattering potential discussed 
in the text. From bottom-to-top, the curves increasingly favor small 
angle scattering and correspond to A z = A^ = <x>, 1.0,0.8,0.6,0.4, 
respectively. Recall that the scale, in n, over which A„ approaches 
/loo corresponds to the real-space length scale of the scattering poten- 
tial within the quasi-2D layer, measured in lattice constants (i.e. it 
corresponds to AT 1 ). The overall magnitude of the increase A^/Ao 
depends on the length scale perpendicular to the layers (i.e. on A7 1 ). 



ity with width in q z c and cp equal to V2A Z and V2A0 respec- 
tively. That is, we take the scattering probability P(\q z \, \<p\) °c 
exp[-(<? z c/2A,) 2 ] exp[-(0/2A ) 2 ] so that 



A„ — A, 



3 (l-e-'^- A 



(38) 



and /lo = /t M (l - e -). This form is used for simplicity, but 
captures the qualitative characteristics of a scattering potential 
with two spatial length scales, one within the layers (K^kfY 1 
and one perpendicular to the layers cAr 1 , for which the scat- 
tering probability is peaked at zero-momentum-transfer. (This 
is illustrated in the Appendix, in which a familiar model for a 
random impurity potential is considered and the collision pa- 
rameters determined to have the same qualitative properties as 
the Gaussian model used in this section.) A plot of the colli- 
sion parameter A n versus its index n is shown, for parameter 
values A-, A^ used below, in Fig. [2] 



1. Small-angle scattering in a 2D system 

In Fig. [3] we illustrate the difference between the behavior 
of AMRO predicted by Eq. [16] and that obtained within the 
relaxation time approximation Eq. [27] We choose parame- 
ters that describe a situation in which small-angle scattering 
is dominant: A- = A^ = 0.3 so that A^/Aq a 10 in order to 
best illustrate qualitative behavior. Also, we use kjc - 3 and 
£2cMo = 15. The curve clearly shows that the relevant scat- 
tering rate in AMRO makes a gradual transition from Aq to A^, 
with increasing 9b- 
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FIG. 3: (Color online.) Going beyond the relaxation- time relaxation 
has a significant effect on AMRO when small-angle scattering is 
dominant. The interlayer resistivity, obtained from Eq. [16] of a quasi 
2D metal as a function of the angle 6 B between the magnetic field and 
the current is shown as the solid line. Both dashed lines are obtained 
using the relaxation time approximation, i.e. by replacing the full 
collision functional with a single scattering rate equal to either the 
transport relaxation rate A or the total quasiparticle scattering rate 
A„. The parameters have been chosen such that A„/Aq ~ 10, which 
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FIG. 5: (Color online.) High-sensitivity of AMRO to momentum- 
dependent scattering in a quasi- ID system. Main panel: The inter- 
chain resistivity in a quasi- ID system is plotted as a function of mag- 
netic field angle for varying scattering parameters with large-angle 
(intersheet) scattering favored. From bottom to top, the curves are 
for scattering that increasingly favors direct scattering from one quasi 
ID Fermi sheet to the other (without momentum change along the 
sheet). The parameters A,, A^ are the same as in Fig. [4] and we have 
used 2tbC/\>f = 0.3. Inset: The solid curve is the A ; = A^ = 0.4 curve 
from the main panel, plotted over the entire range of 9 B , along with 
the relaxation-time results for A s „ = A d n = Aq (lower dashed curve) and 
A s n = A d n = A m (upper dashed curve). 
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FIG. 4: (Color online.) Suppression of AMRO by the removal of 
large-angle scattering processes. The interlayer resistivity of a quasi- 
2D system is plotted as a function of magnetic field orientation for 
varying scattering parameters with small-angle scattering favored. 
Following the arrow, the curves are for scattering that increasingly 
favors small-angle processes, using the same parameters as in Fig. 
|2] The current relaxation rate Aq is the same for all curves while the 
total quasiparticle scattering rate /!„ increases by a factor of seven 
from top to bottom. 



In Fig. |4] we plot the resistivity for various values of 
the interlayer and intralayer parameters Aj. and A$, taking 
kfc - 3 and Q.c/Aq = 15 for all the curves, and use A, = 
A^ = oo,l,0.8,0.6,0.4 from the bottom curve to the top 
along the arrow shown. Note that the transport relaxation rate 
Ao/Q.c = 1/15 is the same for all curves but the total quasipar- 
ticle scattering rate /loo/Qc increases from 1/15 to nearly 1/2 
with ascending curves. The resistivity at n/2 changes little, as 
expected in the strong-field regime. The n-dependence of the 
collision parameters has the effect of increasingly widening 
and suppressing the AMRO as 9b is increased. 



2. Large-angle scattering in quasi- ID systems 

If scattering by a spin or charge density wave with a finite 
ordering wavevector is important, then the scattering proba- 
bility could be peaked about a particular large angle. For the 
quasi-2D system, the isotropic model used here is not applica- 
ble to this case (such scattering only affects electrons near 'hot 
spots' so the system is necessarily anisotropic). However, we 
may look at the case of a quasi-lD system with a scattering 
mechanism that strongly favors intersheet scattering. 

If we ignore scattering within a single Fermi surface sheet 
and take the intersheet scattering probability to have the Gaus- 



sian form above then we find A 

yl _ 1 /1 , „-« 2 a5-a? 



-« 2 M-a 2 ^ 



A a n =,Ul+e- 



A x {\ - e-"- fl ^) and 
). While the former may become small 



for small n, the latter never differ significantly from Aoo- This 
is because direct intersheet scattering, with no change in mo- 
mentum along the sheet, does not change the sum of the elec- 
tron densities of the two sheets but is effective at relaxing any 
difference between them. 

We show, in the main panel of Fig. |5]the resulting magne- 
toresistance for the same scattering parameters as used in Fig. 
|4] Here, we take the parameter 2tbc/(Hvf) = 0.3, which re- 
flects the fact that the interchain hopping parameter % is small 
compared to the intrachain Fermi energy hv/kf. This small 
parameter, which occurs in the argument of the Bessel func- 
tion, restricts AMRO to large 9b- For this reason, only 9b 
between 60° and 90° are plotted in the main panel. The factor 
oJc/'ki — 5 at 9b — as for the quasi-2D case above. 

The resistivity at 9b = 90° increases significantly as the 
parameters A z , A^ are decreased. This dependence occurs be- 
cause the small factor 2tbc/(fiVf) ensures that the second term 
within the square root of Eq. [5T| is not dominant as it is for 
the quasi-2D case. Because of this effect, which magnifies the 
field-angle dependence at large 9b, AMRO are more promi- 
nent for curves with small A,, A^ (i.e. for curves in which 
direct intersheet scattering, without momentum transfer along 
the sheet, is dominant). In contrast, for the A z = A^ = oo 
curve, AMRO are barely perceptible even though the system 
is in the strong-field limit at 9b - 0. 

In the inset, the A c = A^s = 0.4 curve is replotted over the 
full 9b range as the solid curve. The lower and upper dashed 
curves show the result obtained using the relaxation-time ap- 
proximation with tZ - Ao and t b 1 = /t M , respectively. 



0.25 



0.15 




0.05 



FIG. 6: The frequency-width of resonance peaks in inteiiayer re- 
sistivity are determined by corresponding Fourier component of 
the angle-dependent scattering probability. Plotted is the interlayer 
conductivity tr(to), from Eq. [16] at a frequency to = 24/fo ver- 
sus magnetic field strength B. The conductivity has peaks occur- 
ring when the to is an integer multiple of the cyclotron frequency 
a>c = eB/(mcos9 B ). We show the plot for field angle 9 B = 30° 
and have defined a B^ by a> = eB^/m. The peaks at to = u> c ,2a> c 
and 3a>c are visible. Each peak is approximately Lorentzian, with a 
width related to the corresponding collision parameter: i.e., the re- 

hi 



3. Finite-frequency conductivity 

To illustrate the qualitative frequency dependence of the 
conductivity Eq. [16] (or, equivalently, the magnetic field- 
strength dependence of the conductivity at a given finite fre- 
quency) we show, in Fig. [6] a plot of rr(a>) versus field B at a 
very large frequency to = 24Aq and at an angle 9b = 30°. As 
above, we define a field B^ by to — eB^/m*. The plot is for the 
2D system with Gaussian collision parameters: A z = A^ = 1, 
which corresponds to relatively weak 71-dependence in A„ and 
gives Ai/Aq = 1.36 and A2/A0 = 1.55. 

The plot of cr(co) has peaks located whenever the frequency 
a> is an integer multiple of the cyclotron frequency a>c - 
eB cos 6b /m*, i.e. peaks at B cos 9/B^ = 1, 1/2, 1/3, ..etc. The 
height of these peaks drops off rapidly when 9b is small since 
the amplitude of the the pth peak is proportional to a factor 
J 2 p {kfC tan 9b). So, when the argument of the Bessel function 
is of order unity, only the first few peaks are present. The 
peaks are approximately Lorentzian with the width of the peak 
occurring at a> — pa>c being A p /p cos 9b- 

To see the evolution of o-(lo) with angle we show, in Fig. 
Qplots of cr(w) over the same field range at two frequencies: 
a> - 6/I0 in the left panel and u> - 12loq in the right panel, for 
various angles: 9b = 5, 15,30°. The collision parameters are 
the same as in Fig. [6] 

As 9b is increased, the position of the peak in the conduc- 
tivity B - B^l cos 9b moves to higher fields. The peak height 





FIG. 7: (Color online.) The inteiiayer conductivity cr(u>) as a func- 
tion of field B. The frequency in the left panel is half that in the 
right, and the curves are for different field angles 9 B : from bottom 
to top 9 B = 5°, 15°, 30°. The marked peaks in the left panel oc- 
cur where a> = <±>c with u>c = (eB/m) cos 8 B , and those in the right 
panel occur when to = 2to c . In the left inset a zoom of the 6 B = 5° 
curve near to = a>c is compared to a Lorentzian (dashed curve) with 
width Ail cos 6 B . In the right inset a zoom of the 8 B = 30° is shown 
with a Lorentzian of width A 2 /(2 cos 9 B ). The figure illustrates how 
individual collision parameters can be extracted from the inteiiayer 
conductivity at finite frequency. 
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also increases. We noted above that, as long as 9b is not too 
large, the peak in the conductivity where pcj c - u> will be 
a Lorentzian with width A p /(p cos 9b)- However once 9b be- 
comes big enough that J\{kjc tan 9b) I Jo(k)fc tan 9b) « 1, the 
peak will no longer be a single Lorentzian and its width will 
be determined by an average over more than one A p . 

In the left panel, a prominent peak occurs for 9 B as small as 
5°. The inset of this panel shows a zoom view of the 9b - 5° 
conductivity curve near its peak compared to a Lorentzian fit 
with width equal to A\ / cos 5°. Clearly, the curve can be well- 
described by a single Lorentzian and the value of A\ could, in 
principle, be extracted from this type of analysis. In the right 
panel, the peaks are much weaker at small 9b because they are 
associated with the Jiij) term in Eq. Q2] which goes to zero 
rapidly as 9b decreases. This means that curves for larger 9b 
must be considered. For 9b - 30°, a peak is evident and can be 
well fit by a single Lorentzian with width /l2/(2cos30°). At 
this angle there will be slight mixing of the At, and A\ terms 
since [Jx(y)IMy)] 2 « 0.3 at 9 B = 30°. 



V. CONCLUSIONS 

In this article, we calculated AMRO for quasi 2D and ID 
metals using an arbitrary elastic collision integral. The mean- 
ing of the effective relaxation rate in AMRO changes depend- 
ing on the angle between the magnetic field and the current. 
When the field is perpendicular to the layers (and parallel to 
the current), the effective scattering rate is the transport relax- 
ation rate Aq - t^ 1 while for fields along the layer (perpendic- 
ular to the current) the relevant quantity is the total quasiparti- 
cle scattering rate A,*, - r . For intermediate angles, scatter- 
ing is dependent on the in-layer momentum dependence of the 
scattering cross-section, and the associated scattering rate is a 
weighted average of the intralayer Fourier components of the 
scattering probability. We have described methods by which 
these parameters may be extracted, thus allowing detailed in- 
formation about the momentum-dependent scattering rate to 
be obtained. 

It is apparent that the simplicity of the models considered 
here make difficult direct comparisons with data on most sys- 
tems of interest. For example, we considered a 2D system 
that is isotropic in the metallic layers, whereas it is the strong 
anisotropy in the plane that is the focus of interest in many 
2D metals under current investigation. Nonetheless, the qual- 
itative results presented here should be useful for experimen- 
talists in interpreting their data and, in particular, in com- 
paring AMRO scattering rates with independent determina- 
tions. In future theoretical work the combined effects of strong 
anisotropy in the layer and strong momentum-dependent scat- 
tering probabilities should certainly be investigated. 
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VII. APPENDIX: SCATTERING BY A RANDOM STATIC 
POTENTIAL 

To get a clearer picture of the correspondence between the 
collision parameters and the spatial length scales of the scat- 
tering potential, we consider the example of a scattering by a 
dilute, random distribution of impurities^. 

We write the total potential, due to all impurities, at a given 
position f as L^(f) = £ R . u(r - R,) where the sum is over the 
Ni mp impurities in the sample. (The position f is a single- 
electron operator whereas each impurity position R, is a con- 
stant vector since all motion of the heavy impurity is ignored.) 
Using Fermi's Golden rule, the probability Wk,k' that an elec- 
tron is scattered from k to k' is 



Wiu- = ^5(e k - *)l < k|f/(f)|k' > | 2 . 
n 



(39) 



We average this expression over all impurity configurations 
by integrating independently each impurity position R, over 
the sample volume Q. and keep terms to lowest order in the 
impurity density n unp - Ni,„ p /Q.. This gives 



W kM , 



Inn;, 



tiQ. 



-<5(ek-ek')"qM-q 



(40) 



where q = k' - k and u q — J dr exp(-/q • r)u(r) is the Fourier 
transform of the single impurity potential. Comparing this to 
Eq. [3] we obtain the scattering probability P(k, k') as 



P(k, k') 



H(2n) 



2"q M 



(41) 



The probability is related to the collision parameters A mn ac- 
cording to Eq. [T3] If both k and k' lie on the cylindrical 
Fermi surface then the probability depends only on \q z \ and \tp\, 
where (f> is the angle between k and k', and may be expanded 
as P{q z , <f>) - Yimn Pnm exp(imq z c) exp{in(f>). The collision pa- 
rameters Am, are 






with 



Pnm J dq z d<j> exp(-imq z c) e.xp{-inqb)\u(q z , cp)\ 2 

^0 ~ 



fdq z d<f>\u(q z ,(p)\ 2 



(42) 



(43) 



and are thus determined by the Fourier transform of the single- 
impurity potential. 

A simple model of the potential due to a single impurity at 
the origin in a quasi-2D isotropic system is 

u(x,y,z) = u exp(-A z \z\/c)exp(-2A^kf yjx 2 +y 2 ) (44) 

where uq is a constant. The range of the potential within a 
layer is of order (A^fc/) -1 , where kf is the Fermi wavevector 
and A^ is a dimensionless scale. The interlayer range is cA~ l , 
where c is the interlayer lattice constant and A z is dimension- 
less. We have 46 

Arf, 

-n • (45) 

+ sin 2 fp/2 



_ 2 A, 

u(q z ,<f>) = 7TUo(ck f )— — - — -r — 
1 A? + (^c) 2 [ A 2 



11 



If the range of the potential is considerably longer than the 
atomic spacing (so A^A. << 1) then the scattering probability 
is peaked at zero-momentum transfer and begins to decrease 
once the intralayer momentum transfer <p exceeds A^ or the 
interlay er momentum transfer q z c exceeds A e . Considering 
Eq. @2] this implies that A„ m will be much smaller than Aa, if 
m « At 1 and n « A, 1 (the arguments of the exponentials 
are always small and the two terms in the equation nearly can- 
cel in this limit). However, once m,n approach A7 1 and AT 1 , 
the second term in Eq. [42] will start to drop off and thus the 
collision parameters A mn will begin to approach /l M . 

It thus becomes clear that the scale in to, n over which A„,„ 
approaches A*, gives the spatial range of the scattering poten- 
tial perpendicular and parallel to the layers respectively. We 
found out above that only to = 1 terms enter the expression for 
the interlayer conductivity, hence the range of the interlayer 
potential cannot be obtained in this manner. The range of the 



potential within the layers can be obtained, since multiple A\ n 
parameters may be extracted from the interlayer conductivity, 
as discussed above. 

Also, one can at least tell whether the interlayer potential 
extends over a range of significantly more than one lattice con- 
stant by comparing the magnitude of the transport relaxation 
rate A\§ to that of the total quasiparticle scattering rate A\oa, 
both of which can be observed in the interlayer conductivity. 
For, Aio can be significantly smaller than A loo only if the factor 
in the argument of the exponential in Eq. [42j q z c for to = 1, is 
much smaller than unity whenever the scattering probability 
is non-zero (this is the requirement that the two terms nearly 
cancel). This implies that the range of the potential in real 
space A7 1 is much larger than unity. Thus the magnitude of 
the difference between A\q and A\<x, gives a clue as to the range 
of the interlayer scattering potential. 
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